Molecular wires: tuning of electron transport 
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Abstract 

Electron transport characteristics through molecular wires are studied by using the Green's function 
formalism. Parametric calculations are performed based on the tight-binding model to investigate the 
transport properties through the wires. The transport characteristics are significantly influenced by (a) 
the interference effects, (b) chemical substituent group, (c) molecule-to-electrode coupling strength and 
(d) the gate voltage, and, here we focus our results in these aspects. In this article we also discuss the noise 
power of current fluctuations. The noise power gives key information about the electron correlation which 
is obtained by calculating the Fano factor (F) and the complete knowledge of the current fluctuations is 
very essential to fabricate efficient molecular devices. 
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1 Introduction 

The present miniaturization has provided us to ex- 
plore electron transport on the scale of a single 
molecule. The advantage of molecules as an elec- 
tronic building blocks is that the molecules can be 
engineered to have some built in functionality, act- 
ing, for examples, as light sensitive switches, gates, 
or transport elements. Furthermore, molecules are 
very small and hence molecular systems could pro- 
vide a way to down the scale of electronic devices 
even further, especially in self-assembly can be used 
to fabricate nanoscale circuits. Therefore, the field 
of molecular electronics is receiving increasing at- 
tention from fundamental scientists and industry 
alike. In 1974, Aviram and Ratner p] first de- 
scribed electron transport through single molecule 
electronic device-the molecular rectifier. Later, the 
developments of nanoscience and technologies have 
provided several possible route for the construction 
and characterization of single- molecule devices [21 El 
HI 02 El [7] . From experimental developments, the- 
ory can give a better insight in understanding the 
new mechanism of conductance through molecules 
placed between two non-superconducting electrodes 
with few nanometer separation, yet the complete 
knowledge of the conduction mechanism in this 
scale in not well understood even today. There 
are several important factors that control the elec- 
tron transport in such molecular devices. First one, 
of course, is the quantization of energy levels as- 
sociated with the identity of the molecule itself. 
Second one is the quantum interference effects of 
electron waves M EB ED EH EH E] associated 
with the geometry that the molecule adopts within 
the junction. Third are the different parameters of 
the Hamiltonian that describe the molecular sys- 
tem, the electronic structure of the molecule and 
the molecular coupling with the side attached elec- 
trodes. The study of structure-conductance rela- 
tionships is very important for fabrication of effi- 
cient molecular devices with specific properties and 
in a very recent work Ernzcrhof et al. [15j have pre- 
sented a general design principle and performed sev- 
eral model calculations to demonstrate the concept. 
The knowledge of current fluctuations (of thermal 
or quantum origin) in molecular devices is also an 
important issue. Blanter et al. [16] have studied 
elaborately how the lowest possible noise power of 
the current fluctuations can be determined in a two- 
terminal conductor. The steady state current fluc- 
tuations, the so-called shot noise, is a consequence 
of the quantization of charge and it can be used to 



obtain information on a system which is not avail- 
able directly through conductance measurements. 
The noise power of the current fluctuations provides 
an additional important information about the elec- 
tron correlation by calculating the Fano factor (F) 
which directly informs us whether the magnitude of 
the shot noise achieves the Poisson limit (F = 1) or 
the sub-Poisson (F < 1) limit. 

Several ab initio methods are used for the calcu- 
lation of conductance [17l H3 [HI [20l EH [22] through 
a molecular bridge system. At the same time the 
tight-binding model has been extensively studied 
in the literature and it has also been extended to 
DFT transport calculations [23l [24]. The study 
of static density functional theory (DFT) [25l [26] 
within the local-density approximation (LDA) to in- 
vestigate the electronic transport through nanoscale 
conductors, like atomic-scale point contacts, has 
met with nice success. But when this similar theory 
applies to molecular junctions, theoretical conduc- 
tances achieve larger values compared to the ex- 
perimental predictions and these quantitative dis- 
crepancies need extensive study in this particular 
field. In a recent work, Sai et al. [27] have predicted 
a correction to the conductance using the time- 
dependent current-density functional theory since 
the dynamical effects give significant contribution 
in the electron transport, and illustrated some im- 
portant results with specific examples. Quite sim- 
ilar dynamical effects have also been reported in 
some other recent papers ES ES], where authors 
have abandoned the infinite reservoirs, as originally 
introduced by Landaucr, and considered two large 
but finite oppositely charged electrodes connected 
by a nanojunction. Our aim of the present arti- 
cle is to reproduce an analytic approach based on 
the tight-binding model to characterize the elec- 
tronic transport properties through some benzene 
molecules (Fig. [T]) and focus our attention on the 
effects of (a) the quantum interference (b) chemi- 
cal substitucnt group (c) the molecule-to-electrode 
coupling strength and (d) the gate voltage in such 
transport. Here we utilize a simple parametric ap- 
proach [201 El EH 133 131 (33 (Ml EZ1 EM W\ for these 

calculations. The model calculations are motivated 
by the fact that the ab initio theories are computa- 
tionally too expensive, while, the model calculations 
by using the tight-binding formulation are computa- 
tionally very cheap and also provide a worth insight 
to the problem. In our present study, attention is 
drawn on the qualitative behavior of the physical 
quantities rather than the quantitative ones. Not 
only that, the ab initio theories do not give any 
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new qualitative behavior for this particular study 
in which we concentrate ourselves. 

The paper is organized as follow. In Section 2, 
we describe very briefly about the methodology for 
the calculation of the transmission probability (T) , 
current (/) and the noise power of current fluctua- 
tions (S) through a molecule sandwiched between 
two metallic electrodes by using the Green's func- 
tion technique. Section 3 provides the behavior of 
the conductance as a function of the injecting elec- 
tron energy, the current and the noise power of its 
fluctuations as a function of the applied bias volt- 
age for the different molecular wires. In this con- 
text we also discuss the effect of the gate voltage 
on the electron transport through such molecular 
wires. Finally, we summarize our results in Section 
4. 

2 The molecular system and 
the theoretical formulation 

This section describes the models for the molecular 
systems and the methodology for the calculation of 
the transmission probability (T), conductance (g), 
current (/) and the noise power of its fluctuations 
(S) through a molecule (schematically represented 
as in Fig.[T]), sandwiched between the two metallic 
electrodes, by using the Green's function technique. 

Figure [T] shows the schematic representation of 
three different molecular wires in which we con- 
centrate our study. In each wire, the molecule is 
contacted to the two electrodes (source and drain) 
and in addition to that each arm of the molecular 
ring is attached to a voltage gate. The molecular 
system is described by the tight-binding Hamilto- 
nian. Assuming the gate voltages V a and Vb affect 
only on one atom in each arm of the molecular ring 
(atoms a and &), we can write the Hamiltonian for 
the molecule within the non-interacting picture like, 

H M = (eio + V a S ia + VbS ib ) c\ci 

i 

<ij> 

where e^o's are the site energies and t is the nearest- 
neighbor hopping strength. 

At low voltage and low temperature, the conduc- 
tance g of the molecule is given by the Landauer 
conductance formula HDI, 



.9 = ^r T 



where the transmission probability T is written in 
this form BUI. 



T = Tr[T s G M T D G M ] 



(3) 



where G r M {G%[) is the retarded (advanced) Green's 
function of the molecule, and, T$ (Td) describes 
its coupling to the source (drain). The effective 
Green's function of the molecule is expressed as, 



Gm — (E — Hm — Eg — S7 



(4) 



where E is the energy of the injecting electron 
and Hm is the Hamiltonian of the molecule as de- 
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Figure 1: (Color online). Schematic view of 
three different molecular wires where the benzene 
molecules are attached to the two electrodes (source 
and drain) through thiol (S-H) groups in the 
chemisorption technique. The two gate voltages V a 
and Vb are variable. 

scribed in Eq. |T]). In Eq. (@}, £5 and Y, D cor- 
respond to the self-energies due to coupling of the 
molecule to the two electrodes. All the information 
about the molccule-to-clcctrodc coupling are in- 
cluded into these two self-energies and are described 
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through the use of Newns- Anderson chemisorption 
theory [30| [32] . 

The current passing across the molecule can be 
considered as a single electron scattering process 
between the two reservoirs of charge carriers. The 
current-voltage relationship can be obtained from 
the expression [40] . 



I(y) = ^K J (fs-f D )T(E)dE 



where the Fermi distribution function fs(D) = 
/ (E — ns(D)) with the electrochemical potentials 
Ms(-D) = Ep ±eV/2. For the sake of simplicity, here 
we assume that the entire voltage is dropped across 
the molecule-electrode interfaces and this assump- 
tion does not greatly affect the qualitative aspects 
of the I-V characteristics. This assumption is based 
on the fact that the electric field inside the molecule, 
especially for short molecules, seems to have a min- 
imal effect on the conductance-voltage characteris- 
tics. On the other hand, for quite longer molecules 
and high bias voltage, the electric field inside the 
molecule may play a more significant role depend- 
ing on the internal structure of the molecule [41 j . 
yet the effect is too small. 

The noise power of the current fluctuations is cal- 
culated from the following expression [16] , 

oo 

S = ^ J [T(E){f s (l-fs) + f D (l-f D )} 

— OO 

+ T(E) {l-T(E)}(f s -f D ) 2 ] dE (6) 

where the first two terms of this equation corre- 
spond to the equilibrium noise contribution and the 
last term gives the non-equilibrium or shot noise 
contribution to the power spectrum. By calculat- 
ing the noise power of the current fluctuations we 
can evaluate the Fano factor F, which is essential to 
predict whether the shot noise lies in the Poisson or 
the sub-Poisson regime, through the relation [16j . 



F 



S 

2el 



(7) 



For F = 1, the shot noise achieves the Poisson 
limit where no electron correlation exists between 
the charge carriers. On the other hand, for F < 1, 
the shot noise reaches the sub-Poisson limit and it 
provides the information about the electron corre- 
lation among the charge carriers. 

In this article, we present our results at much low 
temperature (5 K), but all the essential features of 



electron transport are also invariant up to some fi- 
nite temperature (~ 300 K). The reason for such an 
assumption is that the broadening of the molecular 
energy levels due to the coupling of the molecule to 
the electrodes is much larger than that of the ther- 
mal broadening. For simplicity, we take the unit 
c=e = /i = lin our present calculations. 



(5) 3 Results and discussion 



Here we describe the electron transport character- 
istics through the three different short molecular 
wires, schematically represented in Fig. [1] where 
the molecules are attached to the two electrodes, 
namely source and drain, respectively. In actual 
experimental set up these two electrodes are con- 
structed from gold and attached to the molecule via 
thiol (S-H bond, i.e., sulfur-hydrogen bond) groups 
in the chemisorption technique, where hydrogen (H) 
atoms remove and sulfur (S) atoms reside. For the 
bridge system given in Fig. HJa), the chemical sub- 
stituent free benzene molecule is attached to the 
two electrodes symmetrically, i.e., the upper and 
the lower arms of the molecular ring have equal 
length. On the other hand, for the bridge given 
in Fig. QJb), a chemical substituent group {CH 3 ) 
is added in one arm of the molecular ring keeping 
the electrodes at the same positions as the previous 
one. Here the symmetry is broken by introducing 
this chemical substituent group which effects the in- 
terference conditions. Finally, in the rest molecular 
bridge, given in Fig.QJc), the electrodes are coupled 
asymmetrically (i.e., the upper and the lower arms 
of the molecular rings have different lengths) with 
the chemical substituent free benzene molecule in 
order to reveal the interference effects on electron 
transport much more clearly. All these three molec- 
ular bridges show various interesting features on the 
electron transport and we will see that the two gate 
voltages, V a and Vb, have significant effect on such 
transport characteristics. 

We will study the behavior of the electron trans- 
port through the molecules in two distinct regimes. 
One is the so-called weak-coupling regime defined 
as T S(D) << £j an d the other one is the strong- 
coupling regime mentioned as t~s(d) ~ t, where ts 
and td are the hopping strengths of the molecule 
to the source and drain, respectively. In our cal- 
culations, the parameters in these two regimes are 
chosen as ts = Tp, = 0.5, t = 3 (weak-coupling) 
and Ts = Tp, = 2.5, t = 3 (strong-coupling). The 
hopping integral in the two electrodes is taken as 
v = 8. 
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Figure 2: (Color online). Conductance g as a function of the injecting electron energy E in the limit of 
weak molecular coupling, where the first, second and third rows correspond to the results for the wires 
given in Figs, [lja), (b) and (c), respectively. The first, second and third columns represent the results 
considering the gate voltages as: V a = 0, Vb = 0; V a = 3, Vb = and V a — 6, Vb — 0, respectively. 



In Fig. [2] we plot the conductance g as a function 
of the injecting electron energy E for the three dif- 
ferent molecular bridges in the limit of weak molec- 
ular coupling, where the first, second and third rows 
correspond to the results for the bridges given in 
Figs.[]Ja), (b) and (c), respectively. The first col- 
umn gives the results for the bridges in the absence 
of any gate voltage (V a = 0, Vb = 0), while the sec- 
ond and third columns represent the results for the 
gate voltages V a = 3, Vb = and V a = 6, Vb = 0, re- 
spectively. From all these curves (of Fig. [2]) it is ob- 
served that the conductance shows very sharp res- 
onant peaks for some particular energies, while for 
all other energies the conductance vanishes. At the 
resonant peaks where the conductance g reaches the 
value 2, the transmission probability T goes to unity 
since we get the relation g = IT from the Landaucr 
conductance formula (see Eq.([2]) with e = h = 1 in 
our present treatment). The resonant peaks in the 
conductance spectra are associated with the energy 
eigenvalues of the single benzene molecules. There- 
fore, the conductance spectrum manifests itself the 
electronic structure of the molecule. From the re- 



sults it is noted that some of the resonant peaks do 
not achieve the value 2 anymore and also get much 
reduced values. This behavior can be explained in 
the following way. The electrons are carried from 
the source to the drain through the molecules and 
the electron waves propagating along the two arms 
of the molecular ring may suffer a relative phase 
shift between themselves. Accordingly, there might 
be constructive or destructive interference due to 
superposition of the electronic wave functions along 
the various pathways. Therefore, the probability 
amplitude of the electron across the molecule be- 
comes either large or small. The ant i- resonances in 
the transmission (conductance) spectra are due to 
the exact cancellation of the transmittances along 
the two paths. The other key observation is that the 
positions of the resonant peaks in the conductance 
spectra get modified with the application of the gate 
voltages. Thus one can get the on/off state of the 
molecular bridge for any fixed energy value or ap- 
plied bias voltage by tuning the external gate volt- 
ages, without changing the structure of the molecule 
itself. This phenomenon is quite significant for 
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Figure 3: (Color online). Conductance g as a function of the injecting electron energy E in the limit of 
strong molecular coupling, where the first, second and third rows correspond to the results for the wires 
given in Figs. HJa), (b) and (c), respectively. The first, second and third columns represent the results 
considering the gate voltages as: V a = 0, Vb = 0; V a = 3, Vb = and V a = 6, V& = 0, respectively. 



fabrication of efficient molecular gates or switches. 
Thus the electron transmission is strongly affected 
by the quantum interference and can be controlled 
by the molecule-to-electrode interface structure as 
well as the external gate voltages. 

Now we describe the behavior of the conductance 
j as a function of the injecting electron energy E 
for these wires in the limit of strong molecular cou- 
pling. The results are shown in Fig. [31 where the 
figures in the different rows and columns correspond 
to the same molecular bridges as in Fig. [2j In 
this strong molecular coupling limit, all the reso- 
nant peaks get broadened substantially compared 
to the weak-coupling case. The strong molecular 
coupling broadens the molecular energy levels and 
accordingly, the resonant peaks become wider (con- 
tribution comes from the imaginary parts of the 
two self-energies, £5 and £d [40]). In the strong- 
coupling limit, the molecular wires conduct elec- 
trons for the wide range of energies compared to the 
weak-coupling case where fine tuning in the energy 
scale is necessary to get the electron conduction 
across the molecules. Thus we can emphasize that 



the molecular coupling strength has a strong influ- 
ence on the electron transport. The other impor- 
tant feature appears from the conductance spectra 
is the existence of the conductance (transmittance) 
zeros. Such anti-resonances are specific to the in- 
tcrferometric nature of the scattering states and do 
not occur in the usual one-dimensional scattering 
problems involving potential barriers. These anti- 
resonant states also appear for such wires in the 
limit of weak molecular coupling but cannot follow 
clearly from the curves (plotted in Fig. [2|) since the 
widths of all the resonances are extremely small. It 
is observed that the positions of the anti-resonant 
states on the energy scale are independent of the 
molecule-to-electrode coupling strength. Since the 
width of these states arc very small, they do not give 
any significant contribution to the current-voltage 
characteristics. However, the variation of the in- 
terference conditions have strong influence on the 
magnitude of the current flowing through the bridge 
systems. From the results we can clearly observe 
the shift of the resonant peaks with the gate volt- 
ages (second and third columns of Fig. [3j. 
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Figure 4: (Color online). Current I and the noise power of its fluctuations S (dotted curve) as a function 
of the applied bias voltage V in the limit of weak molecular coupling, where the red, green and blue 
curves correspond to the results for the wires given in Figs. QJa), (b) and (c), respectively. The first, 
second and third columns represent the results considering the gate voltages as: V a = 0, Vb = 0; V a = 3, 
Vb = and V a = 6, Vb = 0, respectively. 



The scenario of electron transfer through the 
molecular junction can be understood much more 
clearly from the current-voltage (I-V) characteris- 
tics. Here we shall describe the behavior of the cur- 
rent I and the noise power of current fluctuations S 
as a function of the applied bias voltage V for these 
molecular wires, where both the current and the 
noise power are evaluated by the integration pro- 
cedure of the transmission function T (see Eqs.© 
and j6])). Figure [4] displays the current I and the 
noise power of its fluctuations S (dotted curve) of 
the molecular wires in the limit of weak-coupling, 
where the red, green and blue curves correspond to 
the results for the wires given in Figs.QJa), (b) and 
(c), respectively. The first column represents the 
results for the wires in the absence of any gate volt- 
age (V a = 0, Vb = 0), while, the second and third 
columns denote the results for the gate voltages 
V a = 3, Vb = and V a = 6, Vb = 0, respectively. 
Several interesting results appear from these curves 
which we will now discuss one by one. (I) The cur- 
rent (I) shows staircase-like structure with sharp 



steps as a function of the applied bias voltage (V). 
This is due to the sharp resonances those appear 
in the conductance spectra (see Fig. [5]) in the weak 
molecular coupling limit, since the current is evalu- 
ated from the integration procedure of the transmis- 
sion function T. With the increase of the applied 
bias voltage, the electrochemical potentials on the 
electrodes are gradually shifted and eventually cross 
one of the molecular energy levels. Accordingly, a 
current channel is opened up and a jump in the I-V 
curve appears. (II) The current amplitude through 
the molecular bridges strongly depends on the ge- 
ometry of the bridge, which is clearly observed from 
the curves plotted in the first column of Fig. [4] For 
the same bias voltage V the current amplitude for 
the symmetric bridge (Fig. Ufa)) is larger compared 
to the other two asymmetric bridges (Figs.[TJb) and 
(c)). This is due to the quantum interference ef- 
fects of the electron waves traversing through the 
different arms of the molecular ring. The other im- 
portant observation is that for a particular bridge 
the current amplitude can be controlled very nicely 




Figure 5: (Color online). Current / and the noise power of its fluctuations S (dotted curve) as a function 
of the applied bias voltage V in the limit of strong molecular coupling, where the red, green and blue 
curves correspond to the results for the wires given in Figs. QJa), (b) and (c), respectively. The first, 
second and third columns represent the results considering the gate voltages as: V a = 0, Vb = 0; V a = 3, 
Vb = and V a = 6, Vb = 0, respectively. 



by applying the two external gate voltages V a and 
Vb, which is clarified from the results given in the 
second and third columns of Fig. 0] The most sig- 
nificant result is that the threshold bias voltage, 
the voltage where the electron starts conduction 
through the molecule, can be tuned in a control- 
lable way by these two external gate voltages. Thus 
we can tune the current amplitude as well as the 
threshold bias voltage through a molecular bridge, 
externally, by means of these gate voltages. These 
results provide key informations for the fabrication 
of molecular devices. (Ill) In the determination of 
the noise power of the current fluctuations (S) (dot- 
ted curves of Fig. [4]) for these molecular wires it is 
observed that the shot noise goes from the Pois- 
son limit (F = I) to the sub-Poisson limit (F < 1) 
as long as we cross the first step in the current- 
voltage characteristics. This emphasizes that the 
electrons are correlated after the tunneling process 
has occurred. Here the electrons are correlated only 
in the sense that one electron feels the existence of 
the other according to the Pauli exclusion principle, 



since we have neglected all other electron-electron 
interactions in our present treatment. 

Now focus our attention on the current and the 
noise power of its fluctuations for the molecular 
wires in the limit of strong coupling. The results 
are shown in Fig. [5l where the red, green, blue and 
dotted curves correspond to the same meaning as 
in Fig. 3J In this strong coupling limit we also 
get several important features and here summarize 
them. (I) The current varies almost continuously 
with the applied bias voltage for all such molecu- 
lar wires. The key point is that the current ampli- 
tudes get enhanced quite significantly compared to 
the weak molecular coupling limit. This behavior 
can be clearly understood by noting the areas un- 
der the curves in the conductance spectra for this 
strong molecular coupling limit (see Fig. [3|). Thus 
for a particular bridge system one can enhance the 
current amplitude by increasing the moleculc-to- 
electrode coupling strength. This is an interesting 
phenomenon in the study of molecular transport. 
(II) Like as in the weak-coupling case, here also 
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the current amplitude decreases due to the break- 
ing of the symmetry of the molecular bridge (see 
the second and third rows of Fig. [5]). The reason 
for such behavior is the same as in our previous de- 
scription. Here also the current amplitudes are con- 
trolled by the external gate voltages. (Ill) Finally, 
in the study of the noise power of the current fluc- 
tuations in this strong coupling case we see that, 
for the molecular bridge given in Fig. [lja) there 
is no such possibility of getting a transition from 
the Poisson limit (F = 1) to the sub-Poisson limit 
(F < 1) since the shot noise already achieves the 
sub-Poisson limit (see the dotted curves of the first 
row in Fig. [5]) , momentarily as we switch on the bias 
voltage. Therefore, for this particular bridge in this 
limit of molecular coupling the electron correlation 
is highly significant. On the other hand, for the 
other two bridges the shot noise makes a transition 
from the Poisson limit to the sub-Poisson limit af- 
ter some critical value of the applied bias voltage 
(see the dotted curves of the second and third rows 
in Fig. [5|). This critical value of the bias voltage 
depends on the geometry of the molecular bridge 
as well as the external gate voltages for the fixed 
molecular coupling strength. 

4 Concluding remarks 

To summarize, we have studied the electron trans- 
port properties through some molecular wires, 
based on the tight-binding model, by using the 
Green's function formalism. The transport proper- 
ties through the wires are significantly affected by 
the several factors like, (a) quantum interference of 
the electron waves traversing through the different 
arms of the molecular ring, (b) chemical substituent 
group, (c) molecular coupling to the electrodes and 
(d) the external gate voltages. All the characteris- 
tic features described here provide several key ideas 
for fabrication of efficient molecular devices. 

For the weak molecular coupling limit the con- 
ductance shows fine resonant peaks (Fig. [2|) for 
some particular energy values, while, for all other 
energies it drops to zero which are in fact the sig- 
nature of the electronic structure of the molecules. 
On the other hand, in the limit of strong molecu- 
lar coupling the width of these resonant peaks get 
broadened substantially (Fig. [3|) . This is due to the 
broadening of the molecular energy levels, where 
the contribution comes from the imaginary parts of 
the self energies Eg and £d [40]. When the molec- 
ular symmetry is broken, more anti-resonant peaks 
appear in the conductance spectra and their posi- 



tions are independent of the molecule-to-electrodc 
coupling strength. 

The scenario of the electron transfer through 
the molecular bridges can be visible much more 
clearly by studying the current- voltage characteris- 
tics. Current shows the staircase-like behavior with 
sharp steps (colored curves in Fig. |4|) in the limit of 
weak molecular coupling, while, it gets a continuous 
variation (colored curves in Fig. [5]) as we increase 
the coupling strength and achieves much larger am- 
plitude. The current amplitude can also be tuned 
by applying the external gate voltages (second and 
third columns of Figs. [4] and [5]). 

Finally, in the determination of the noise power of 
the current fluctuations we have seen that whether 
the shot noise lies in the Poisson limit (F = 1) or 
the sub-Poisson limit (F < 1) strongly depends on 
the geometry of the molecular bridge as well as on 
the molecular coupling strength (dotted curves in 
Figs. H and O. 

Throughout our discussion, we have used sev- 
eral approximations by neglecting the effects of 
the electron-electron interaction, all the inelastic 
scattering processes, the Schottky effect, the static 
Stark effect, etc. More studies are expected to take 
into account all these approximations for our fur- 
ther investigations. 
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